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Abstract. - We demonstrate the possibility of a stable equilibrium multi-lamellar ("onion") 
phase in pure lamellar systems (no excess solvent) due to a sufficiently negative Gaussian 
curvature modulus. The onion phase is stabilized by non-linear elastic moduli coupled to a 
polydisperse size distribution (Apollonian packing) to allow space-filling without appreciable 
elastic distortion. This model is compared to experiments on copolymer-decorated lamellar 
surfactant systems, with reasonable qualitative agreement. 



Introduction. - The existence of vesicles at thermal equilibrium is a long-standing, and 
still controversial, problem [1]. Unilamellar and multi-lamellar vesicles (MLVs) were first 
induced in lyotropic lamellar phases with excess water by adding energy (e.g. shear flow, 
ultrasound, electric field) [2]. In a very elegant work, using membranes generated by a chemical 
reaction, Hoffman and co-workers demonstrated that any of lamellae, unilamellar vesicles, or 
MLVs can be prepared in the same system, depending on the mechanical path chosen [3]. 
However in some special cases [4, 5] equilibrium unilamellar vesicles have been demonstrated. 
These systems are all in the dilute regime (large excess of water). Some experiments have 
suggested that MLVs (or onions) can be stabilized in the semi-dilute regime (excess water) [5]. 
Theories to explain the stability of dilute unilamellar vesicles either describe a competition 
between the entropy of mixing and the curvature energy of the vesicles [5-7] or a symmetry- 
breaking instability leading to a spontaneous curvature [8]. Onions are also predicted to be 
stabilized in the dilute and semi-dilute regimes due to an unstable curvature energy; in these 
cases, a transition towards unilamellar vesicles is avoided by imposing either a core energy [9] 
or a cutoff in the entropy of the Helfrich interactions [10]. Indeed, Simons and Cates [10] have 
discussed the stability of unilamellar vesicles and the transition to onions as the concentration 
increases even when the curvature energy is unfavorable, due to entropic reasons. 

In the concentrated regime (homogeneous lamellar phase with no excess water), very 
monodisperse onions can be prepared under shear [11-13]. The vesicles fill space and remain in 
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the one phase region without expelling solvent, by distorting into space- filling polyhedra [14]. 
The resulting texture is a lattice of disclinations which can be either disordered [11] or or- 
dered [13]. Moreover many experimental systems "spontaneously" exhibit onions in regions 
of the phase diagrams that seem to be continuously linked to the lamellar phase [15]. This 
evolution from a texture of polydisperse vesicles to a texture of focal conies upon tuning a 
parameter, such as a co-surfactant concentration or the volume fraction of solvent, seems to 
be generic. A striking signature of this evolution is the change in the rheology, from a viscous 
gel in the "onion" phase to a more fluid phase in the "lamellar" focal conic phase. Thus, 
whether an onion texture is at equilibrium or not seems yet unclear. In this work we propose 
a model based on a non-quadratic elastic energy, which shows that a pure lamellar phase, or 
one of space-filling deformed monodisperse onions, is unstable with respect to a polydisperse 
space-filling packing (we study an Apollonian packing [16]). We compare these predictions to 
experiments on copolymer-doped lamellar phases, originally introduced by Ligoure et al. [17]. 

Model. - The bending free energy per area FJA for membranes is usually taken, following 
Helfrich, as the simplest quadratic function of the mean and Gaussian curvatures, respectively 
2H = ri 1 +r^ 1 and G — l/(nr2), where n and ri are the two principal radii of curvature [18]. 
This suffices for large curvature radii r » 5, where S is the layer thickness. However, if onions 
exist, because of the high concentration of membrane, high curvatures (small r) are expected 
near the core of the onions. We will thus use the free energy, relative to a flat state, of a 
symmetric bilayer expanded to quartic order in the curvature radii [19]: 

j=2kH 2 + RG+ ± Cl # 4 + \c 2 G 2 + 2c 3 GH 2 , (1) 

where k and R are the conventional bending and Gaussian curvature moduli, and the non- 
linear moduli a can be expected to scale as ksTS 2 . 

Stability of a flat membrane requires 2k > — R > 0. For R > an instability towards 
a phase with G < (e.g. a bicontinuous phase) is possible. However, for R = 2k + R < 
spherical shells are stabilized, with free energy (for radius r) F s h e \\(r) < given by 

FehellW =47TK+^-, (2) 

where c = ci + 4c2 + C3. Positive c stabilizes spherical shells (vesicles) of finite radius. 
Unilamellar vesicles cannot accommodate volume fractions close to space filling, so we consider 
the free energy of an onion of radius R, constructed from k shells of discrete radii rj = jd, (j = 
1, 2, 3, . . . , k) set by the layer spacing d, which we assume remains fixed: 

F onion (R = kd)=J2 ^heiifo - 3d) = - J [Xk + S(-2, k)} (3) 

where A = Rd 2 /c is the balance between the Helfrich and non- linear elastic coefficients and 
S(l,m) = Y^iLi^i 1 )- I 11 th° limit R » d the second term above, due to the fourth order 
coefficients, is equivalent to the core energy introduced by Fournier and Durand [9]. 

For large dilute monodisperse onions the free energy per unit volume scales as F/V ~ 
-Fonion(-R) ~ c(A + d/R)/(d 3 R 2 ) 1 leading for A < to stable onions with size R ~ d/\X\. 
However, in the concentrated regime monodisperse onions cannot fill space without non- 
trivial deformations [14]. Typically, onions formed under shear pack as polyhedra, with the 



( 1 )There is generally, a surface energy due to, e.g., van der Waals attraction, but if we assume an inter-onion 
spacing identical to the equivalent lamellar separation wc need not consider such a term. 
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Fig. 1 - (a) Maximum number of layers per onion, fc max , as a function of A = kd 2 /c. Apollonian 
onions are stable for —1.125 < A < —0.86. (b) Schematic phase diagram in the (k, (f>) plane, k* ~ <j) 2 . 
(c) Schematic packing of the Apollonian distribution (reprinted with permission from [16]). 



deformation of layer j concentrated along the edges of length Vj with curvatures of order 1 /£, 
where the smectic bending penetration length £ = y/ K/B depends on the bending K = n/d 
and compression B moduli of the lamellar phase [20] . Summing over all shells in a deformed 
polyhedron leads to an energy per unit volume that scales as F e i/V ~ c[A + d 2 / £ 2 } / (R£d 3 ) . 
For A + cP/£ 2 < this leads to an unphysical unilamellar vesicle phase with R — d (the lower 
cutoff). However, a polydisperse phase can relax the deformation, leading, for a broad enough 
distribution to a space- filling distribution of spheres. A candidate space-filling distribution is 
the Apollonian packing, for which the number distribution of spheres n(R) ~ (R/d)~ DA , with 
a fractal dimension estimated from simulations to be Da — 3.45 [16,21,22] (see Fig.^). Such 
a packing has been previously suggested to describe focal conic defects in smectic droplets 
[23]. In principal we should calculate this distribution from a free energy that incorporates 
entropy. Although small on the typical scale of the elastic terms, this entropy should widen 
the parameter range in which onions are stable [10]. Existing equations of state for, e.g. 
polydisperse hard spheres, fail for widely disperse distributions near space filling [25], so we 
have been unable to estimate this contribution to any useful accuracy. Entropic effects will 
most likely influence the upper and lower cutoffs of the distribution, and would eventually lead 
(for highly fluctuating systems) to an exponential, rather than power law, size distribution 
away from space filling [7,9]. 

The normalized number distribution of onion sizes Rk obeying -y-X)fc=i x n (Rk)R% = Vm> 
where V m is the volume of lamellar material, is 
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The total energy of the distribution, F tot = Ylk=T n{k)Fonian(k), is 
A7rd 3 F tot _ A S(l — Da, U) + S?°r ^(~ 2 > 3)3 

"-max J 



-Da 



(4) 



(5) 



This energy can then be minimized over the maximum onion size fc max as a function of the 
single parameter A to find the lowest free energy distribution (Fig. ^,). Although spheri- 
cal shells are stable for A < 0, the stabilizing elastic constants Cj assure that a finite and 
negative A < A* ~ —0.86 is required for onion formation. Onions are first stable with 'in- 
finite radius', implying a smooth transition from lamellae to onions. However, the scaling 
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Fig. 2 - Top: Gel phase observed between crossed polarizers, for membrane fractions cj> and copoly- 
mer/surfactant weight ratios a of (a) <j> — 16%, a = 0.2; (b) <f> — 16%, a = 0.4; (c) (f> = 16%, a = 0.8; 
(d) (j> — 22%, a = 0.8. Bottom: Gel phase in contact with water; (e, crossed polarizers) <f> — 16%, 
a = 0.15; (f, Differential Interference Constrast microscopy) <j> — 16%, a — 1.6. Scale bar 50 /im. 



nR 3 ~ R 3 ~ Da ~ R~ 0A5 implies an unbounded total volume. Hence, the size distribution 
must have an upper cutoff, which in practice should be the smallest dimension of the sample 
container. (In fact, the nature of the size distribution for Apollonian packings of spheres 
is still under debate [21,22], and appears to depend crucially on the manner in which the 
packing is constructed.) The maximum onion size then decreases, and for A < —1.125 the 
stable phase is of unilamellar vesicles. At this point the calculation has presumably broken 
down, and either another phase such as micelles intervenes ( 2 ), or even higher order curvature 
terms become important. Because A — kd 2 /c, and the volume fraction <ft is of order 5/d, the 
condition A < A* implies a phase boundary k* ~ <fi 2 , as shown in Fig. ^d, in the k~4> plane. 

Experiments. As shown theoretically in [26], adsorption of amphiphilic copolymers 
to surfactant bilayers is expected to decrease k, which eventually becomes negative. We 
therefore prepare lamellar phases decorated with copolymers [17]. Bilayers of thickness Sq = 
2.8 nm comprising cetylpyridinium chloride (CpCl) and octanol (Oct) (CpCl/Oct = 0.95 
w/w) are diluted in brine ([NaCl] = 0.2M) and decorated with Symperonics F68 (Serva) 
[(EO)76 — (PO)29 — (EO)76, where EO is ethylene oxide and PO is propylene oxide]) as am- 
phiphilic copolymer. The bilayer volume fraction (j) and the copolymer/surfactant weight 
ratio a range from 9 to 23% and to 1.6, respectively. The bilayer volume fraction (j> and 
the copolymer/surfactant weight ratio a range from 9 to 23% and to 1.6, respectively. The 
lamellar structure is preserved upon adding copolymer, with a continuous hardening into a 
"lamellar gel" . The location of this gel in the (a,cj>) phase diagram was studied by Ligoure et 
al. [17]. They observed gel-like behavior over a range of critical copolymer/surfactant weight 



( 2 )A cubic phase of unilamellar vesicles has been experimentally observed in the concentrated regime [24]. 
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Fig. 3 - (a) Schematic phase diagram, adapted from [17], overlaid with the experimental points (x), 
where a and <j) are respectively the copolymer/surfactant weight ratio and bilayer volume fraction. 
Experimental phase boundaries are inferred at (•). The dotted line is a guide for the eye. (b) Small- 
and large-angle static light scattering intensity for (j> = 16% and a — 0.8. A crossover from power law 
to constant intensity is noted at a length scale of order 7 pm. 



ratios that depends on the membrane fraction. 

The marked variation of the mechanical properties of the bulk lamellar phase is accompa- 
nied by modifications of the optical textures. The pictures shown in Fig. ^demonstrate that 
the lamellar gels arc onion phases. Onions are clearly obtained for both neat lamellar phases 
and samples put in contact with solvent. In the latter case individual onions detach from 
the solvent/gel interface (Fig. Eh,f), while in the former case defect textures characteristic of 
the onion phase are observed between crossed polarizers (Fig. |2t--c). The onions are highly 
polydisperse, with an apparent maximum size which seems to decrease with increasing a. 
Different textures and the presence of oily streaks are discerned at higher <fi (Fig.[2Jl)- 

While microscopy indicates that the size distribution is very broad and varies with a and cf>, 
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Fig. 4 - Variation of the elastic plateau modulus Go (top) and maximum radius Rmax (bottom) 
determined from Go, for samples with (a) <f> — 16% and various a; (b) a = 0.8 and various (j>. Solid 
lines are fits to R ~ (a - a c )~ p and R ~ (<j) c - <t>)~ m , with a c = 0.08, <f) c = 24.5%, p = 1.3, m = 3.0. 
(Inset, a) same data as in the main figure in a semi-logarithmic plot. 
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a quantitative determination of the size distribution using this technique is delicate. However, 
assuming an Apollonian distribution we can, for example, estimate the macroscopic elastic 
modulus as a mass (or volume) average of the elastic moduli of onions of different sizes r^. The 
elastic energy stored in a linear deformation of strain /i of an onion of radius R is AirjR 2 ^ 2 , 
where the effective surface tension 7 = |V KB [27], leading to a modulus Gr = 3ttj/R. Thus, 
an estimate of the volume-averaged elastic modulus is Go — y J2i^ n i r i [20]- Averaging for 

/ \ i-D A 

Rmax > d = S /cf>, we find Go ~ % J ~ 3.67 7 cT - 45 i? m ° a f for D A = 3.45. 
The linear elastic plateau is measured in a Couette rheometer along two lines in the phase 
diagram at constant a and <fi (Fig.|3Ji). For fixed 0, the modulus Go increases monotonically 
with increasing a, while for fixed a, the modulus varies non-monotonically with ip (Fig.^J. To 
extract the maximum onion radius R,„ aj . from Go we must estimate the surface tension 7. For 
interlamellar forces dominated by Helfrich entropic undulations, V KB = (d-J) 2 > w ith $ the 
bilayer thickness [17,28]. Interestingly, it has been shown [17] that for a decorated lamellar 
phase the inter-bilayer interactions are still of an effective Helfrich form, with an effective 
bilayer thickness b e g = 5q + 2 x h po i where h po i is the apparent thickness of the polymer 
layer ( 3 ). In the experimental range of a we take h po i = Aa 1 / 3, with A — 3.5 nm, although 
this expression is strictly valid only in the brush regime (a > 0.5). One can thus calculate 
the surface tension by taking S = 5 e ff and hence extract, from the experimental values of 
Go, the variations of R m ax with the two experimental parameters <fi and a. As can be seen 
in Fig. 21 for <fi — 16%, R m ax varies between 1 /xm and 120 /xm, with an apparent divergence 
near a c ~ 0.08; while for a = 0.8, R m ax varies between and 1400 /xm, with an apparent 
divergence near <fi c ~ 24.5% ( 4 ). 

The critical-like behavior of the onion maximum size inferred experimentally is remarkably 
similar, qualitatively, to the prediction of Fig.Q] From the observed critical-like variations in 
Rmax we have identified two points on the boundary between lamellar and onion phases, in 
the (a, 4>) plane (Fig.|^i). Using the relation between k and the parameter A we can estimate 
experimental values for k. The onion/lamellae phase boundary is given by A* = = -0.86 
with c = ksTS^ff. For the two copolymer/surfactant weight ratios a for which the boundary 
has been determined (Fig. ^), we obtain k — — O.lfcsT for a — 0.08 and k = — 0.6/csT for 
a = 0.8. In agreement with theoretical expectations [26], k is found to be negative, of order 
1 ksT, and to decrease with increasing weight ratio. Consistent with the theoretical phase 
diagram (Fig. \\]p) and with previous experimental observations [17], we also find that the 
critical value of k, below which the onion phase is stable, decreases when <f> increases. 

The radii R m ax estimated from the elastic moduli are similar to those observed by mi- 
croscopy (Kig. [2]3,f ) ■ For a — 0.15 the largest observable onion is R ~ 35 /im and the calculated 
Rmax(Go) = 120 /im (Fig.|2^), and for a = 1.6, the largest observable onion is R ~ 5/xm (in 
Fig. Et) an d the calculated R m ax{Go) — 1.4 /im. We can also compare R ma x to static light 
scattering experiments (Fig.yJ. The scattering intensity I(q) scales as a power law for large q, 
I(q) ~ q~ 3 ' 7 , and has a plateau for smaller q, indicating that there is no structure at lengths 
larger than the cross-over length 27r/g ~ 7/im. This provides an estimate for the maximum 
onion size R m ax that is very close to the value obtained from Go (Rmax = 6 fim) . While 
similar power law behavior for high q has previously been observed for monodisperse onions 
obtained by shear, the behavior at large q [29] was very different, with a peak characteristic 

( 3 ) Since the lamellar phase is stabilized by Helfrich interactions, one should in principle incorporate layer 
undulations into an onion entropy [10]. However, we believe that these effects, important for small onions in 
dilute systems, would not significantly change the phase diagram at the higher concentrations studied here. 

( 4 ) We note that, at small <j>, R m ax is a decreasing function of tjt, which may be be due to the nearby lamellar- 
to- vesicle phase boundary [17]. 
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of the onions size observed in that case. The flat intensity obtained here clearly indicates the 
intrinsic polydispersity of the onions. Figures f also show wide size distributions. 

In summary, we have presented experimental evidence for an equilibrium phase of polydis- 
perse multi-lamellar vesicles (onions), and rationalized this by combining a non- linear Hclfrich 
elastic theory with an assumed polydisperse distribution of onions. Future work should incor- 
porate the entropy and calculate the nature of the space-filling distribution. 
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